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ABSTRACT 

We report on a blind survey for extragalactic radio variability that was carried out by comparing 
two epochs of data from the FIRST survey with a third epoch from a new 1.4 GHz survey of SDSS 
Stripe 82. The three epochs are spaced seven years apart and have an overlapping area of 60 deg^. We 
uncover 89 variable sources down to the millijansky level, 75 of which are newly-identified, and we find 
no evidence for transient phenomena. This new sample of variable sources allows us to infer an upper 
limit to the mean characteristic timescale of AGN radio variability of 14 years. We find that only 1% 
of extragalactic sources have fractional variability fvar > 3, while 44% of Galactic sources vary by 
this much. The variable sample contains a larger fraction of quasars than a comparable non-variable 
control sample, though the majority of the variable sources appear to be extended galaxies in the 
optical. This implies that either quasars are not the dominant contributor to the variability of the 
sample, or that the deep optical data allow us to detect the host galaxies of some low-z quasars. We use 
the new, higher resolution data to report on the morphology of the variable sources. Finally, we show 
that the fraction of sources that are variable remains constant or increases at low flux densities. This 
may imply that next generation radio surveys with telescopes like the Australian Square Kilometer 
Array Pathfinder and MeerKAT will see a constant or even increasing fraction of variable sources 
down into the submillijansky regime. 

Key words: galaxies: active quasars: general radio continuum 



1. INTRODUCTION 

Radio sources with variable continuum emission cover 
a broad range of interesting astrophysical phenomena, 
from fiaring stars to accreting supermassive black holes. 
Within our Galaxy, we see radio variability from pulsars, 
magnetars, microquasars, brown dwarfs, cataclysmic 
variable s, and many low mass stars of spectral types K 
and M (jBergeij [2001 . Outside the Milky Way, we see 
exotic phenomena like Gamma-Ray Bursts, radio super- 
novae (Wciler et al. 2002), and blazars. Moreover, as the 
time domain is still relatively unexplored, new classes of 
sources continue to be discovered, su ch as the mysteri- 
ous r otating radio transients (RRAT; iMcLaughlin et all 
[20061) . 

Several upcoming surveys are citing the untapped po- 
tential of the time domain and making variability a prior- 
ity. The Australian Square Kilo meter Array Pathfinder 
(ASKAP; iJohnston et al.l I2008D is one upcoming tele- 
scope with variability as a priority. Under development 
in a remote region of Western Australia, it will consist of 
36 antennas with wide field of view phased array feeds. 
The ASKAP Survey for Variables and Slow Transients 
(VAST; iMurphvet al.ll20131 is being planned with the 
goal of characterizing the radio transient and variable 
sky. Apertif is a project to upgrade the Westerbork Syn- 
thesis Radio Telescope (WSRT), increasing its field of 
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view by a fac tor of 25 an d facilitating the detection of 
variables dOosterloo et al.l[2010) . Transient sources are 
one of the prime science drivers for MeerKAT, which 
is being built in the Northern Cape of South Africa 
and will be the largest and most sensitive radio tele- 
scope in th e Southern Hemisphere upon its completion 
(iBooth fc J onas 2012). The Karl G. Jansky Very Large 
Array (VLA) represents an upgrade to the original VLA 
that is more than 10 times more sensitive and will allow 
the quick detection of varible and t ransient sources. The 
Murchison Widefield Array fMWA: [Lonsdale et al.ll2009D 
will be sensitive to tran sient radio events in the range 80- 
300 MHz, and LOFAR (jHessels et al.ll2009l ) has made the 
search for low-frequency (10-250 MHz) transients one if 
its key science goals. Further in the futu re, the Square 
Kilometer Array fSKA: [van Haarleml2000D is expected to 
contribute significantly to the transient parameter space. 
A number of papers, many using the first Allen Tele- 
scope Array (ATA) results, have recentl y been pub 



lished on searches for transient sources (iBower et al 



20071: ICroft et al.l[20Tol: [Bower et al.l[20Tl IBower fc Saul 



20111 ). Studies of long -term variabi l ity ar e harder to 
find. One such study, bv [Becker et al.l (|2010[) . used three 
epochs of VLA observations to look for variable sources 
in the Galactic Plane. Their 6 cm data covered 23 deg^ 
of the Galactic Plane down to a limiting flux density of 
1 mjy. They found 39 variable sources, and they showed 
that these sources are more highly variable than extra- 
galactic objects. They also concluded that the variable 
fraction increases toward the Galactic center. 

Out of the Galaxy, radio-loud active galactic nuclei 
(AGN) are by far the most common object in radio imag- 
ing surveys with fiux thresholds above ^1 mJy. Radio- 
loud AGN have been observed to vary on timescales 
from less than a day to years in all frequency regimes. 
Some of the radio variability, particularly at low (<1 
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Fig. 1. — Area covered by each of the three epochs. Epoch I and II come from the FIRST survey and are shown as blue dots and red 
squares, respectively. Epoch III is from a recent high-resolution VLA survey of Stripe 82 and is shown by the black dots. Note that the 
axes are not to scale. 



GHz) frequencies, may be extrinsic to t he source. For 
example, i nterstellar scintillation (ISS; iRickettl 119861 : 
iHeescheiifc Rickett 1987 ) and extreme scatte ring events 
(ESE: iFiedler et al.llT987l: iRomani et aIlll987D are inten- 
sity variations thought to be caused by the turbulent, 
ionized interstellar medium of our own galaxy. However, 
ESEs are extremely rare, and ISS produces typ ical in- 
tensi t y var iations of only a few percent (e.g., Love ll et al.l 
120031 12008!). A majority of the observed extragalactic ra- 
dio variability is therefore thought to be intrinsic to the 
source itself. 

A popular theory for the instrinsic variability in 
extragalactic radio sources involves shock waves pro- 
pogating along an adiabatic, conical, relativistic jet 
(jMarscher &: GeaillT9 85). In this model, the amplitude of 
variability should be larger (and the timescales shorter) 
for objects viewed close to the line of sight of the relativis- 
tic jet, such as BL Lac ob j ects and flat spectrum q uasars 
(IBlandford fc R.eesI \WM lUrrv fc Padovanil rT995h . Ex- 
tensive work has therefore be en done on the extreme vari- 
ability of bright blazars (e.g..lHufiiagel fc Bregmanlll992l : 
lAller et al.|[l999l:lCiaramella et al.ll2004D . despite the fact 
that such sources make up only a small fraction of the 
total AGN population. Other studies have targeted 
small gr oups of specifically- selected objects, such as the 
study by lFalcke et al.l (|2001| ) of 30 radio-quiet and radio- 
intermediate AGN, or they have concentrated on corre- 
lating shorter-term variability at multiple frequencies. 



Th ere have been surpri singly few blind, long-term stud- 
ies. iCarilli et al.l ()2003| ) used multi-epoch VLA observa- 
tions to study submillijansky variability i n the Lockman 
Hole, producing just nine variable sources. Ide Vries et al.l 
()2004[) used two VLA epochs of Faint Images of the Ra- 
dio Sky at Twenty centimeters (FIRST) data to search 
for var iabi lity dow n to 2 niJy oy er 104 deg^. More re- 
cently, [Thvagaraja n et al.l (|2'oil[ ) used -55,000 FIRST 
survey snapshots to search for variable and transient ob- 
jects on timescales of minutes to years. This last survey 
produced 1,627 variable and transient sources down to 
millijansky levels, the largest sample to date. 

In this paper, we use a new 1.4 GHz survey of the Sloan 
Digital Sky Survey (SDSS) Southern Equatorial Stripe 
f|Hodgc ct al.ll201l[ ) to further our understanding of long- 
term radio variability in faint extragalactic sources. This 
survey of the Southern Equatorial Stripe, a.k.a. "Stripe 
82" , covers 92 deg^ to a typical rms of 52 /iJy, making 
it the widest survey to reach this depth. These high- 
quality data ove rlap with the FIRST coverage of Stripe 
82 examined by Ide Vries et al.l (|2004l ) for variability, al- 
lowing us to extend their pilot study and add a third 
epoch of data with three times the angular resolution 
(1.8"). 

We begin by describing the radio data in Section [5J 
The details of our sample selection are discussed in 
Section |3l We describe our results in Section JH in- 
cluding sections on source strength dependence (Sec- 
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tion l4.ip . fractional variability (Section 14.21) . characteris- 
tic timescale (Section 14.31) . optical counterparts (Section 
14. 4p , and morphology (Section 14.51) . We discuss the re- 
sults in Section [5l and we end with our Conclusions in 
Section [6l Where applicable we assume the standard 
ACDM cosmolo gy of ffn = 7C I km s~^ Mp c~\ n^ = 0.7, 
and flM = 0.3 (iSoergel et all 120031120071) . 



2. RADIO DATA 

We used three different epochs of radio data to conduct 
our study of radio variability. The first two epochs (1995 
and 2002) come from the FIRST survey (|Becker et al.l 
I1995D . The FIRST survey was conducted with the Very 
Large Array (VLA) at 1.4 GHz and covered >9,000 deg^ 
in the North and South Galactic caps. The survey uti- 
lized the VLA's B-configuration, giving it a resolution 
of 5". The typical rms achieved was 0.15 mJy, and the 
catalog used a 1 mJy detection threshold. 

The Southern Equatorial Stripe was first observed by 
the FIRST survey in 1995. The area was subsequently 
re-observed in 2002 for two reasons: as a quality control 
test of the FIRST survey, and to search for variability. 
The area that was re-observed is smaller than the initial 
area observed, covering 104.3 deg^ and containing 9,086 
radio sources. We will refer to the 1995 and 2002 epochs 
as Epoch I and Epoch II, respectively. 

The third epoch of radio data comes from a re- 
cent high-resolutio n VLA survey of SDSS Stripe 82 
(jHodgeet al.ll20Tl[) . This 1.4 GHz survey covered 92 
deg^ to a typical rms of 52 /xJy, making it the largest 
1.4 GHz survey to reach such a depth. The majority of 
the data were taken in the A-configuration, with sup- 
plemental B-configuration data obtained on every field 
to increase sensitivity to extended structure. The reso- 
lution achieved was 1.8", or roughly three times better 
than FIRST. These data constitute Epoch HI and may 
sometimes be referred to as the 'Stripe 82' radio data. 

The area covered by Epoch HI does not entirely overlap 
the area covered in Epochs I and II. Figure [1] shows the 
area covered by the individual epochs. The SDSS Stripe 
82 data (Epoch HI) avoided the region around RA — 0, 
as planned infrared observations would be compromised 
by zodiacal light in that part of the sky. The total area 
o bserved by all three e pochs is 60 deg^. 

Ide Vries et al.l ()2004[ ) previously searched for variables 
in Stripe 82 using the data from Epochs I and II. In the 
course of their analysis, they determined that the flux 
densities of the Epoch II data required a small correc- 
tion factor consisting of a 90 fj,3y zero-point offset and 
a 1.16% sensitivity correction, well within the estimated 
~5% systematic uncertainty in the flux density scale. In 
the following analysis, we have used the Epoch II data 
with this correction applied. 

3. SAMPLE SELECTION 
3.1. Variable Sources 
This work builds on the work that Ide Vries et al.l 



(|2004[) previously did to search for variability between 
the two FIRST epochs (I and I I). Of the 128 variable 
sources found by Ide Vries et al.l (|2004D between Epochs 
I and II, 58 fall in the area covered by the Epoch HI ob- 
servations. Before searching for new variables with the 
Epoch HI data, we first matched these known variable 



sources to the Epoch HI data to look for systematic ef- 
fects between the epochs. 

Using a 3" matching radius, we matched these sources 
to the Epoch HI data. We chose this matching radius 
as it is approximately half of the FWHM of the beam 
from the two earlier FIRST epochs. We found that 57 of 
the 58 sources were recovered in the Epoch HI catalog, 
with the single source that was not recovered falling just 
below the detection threshold for the Epoch HI catalog. 
We obtained the peak and integrated flux density of this 
source by using the task JMFIT in AIPS. 

Figure[2^ shows the Epoch I/II peak flux density of the 
previously known FIRST variables plotted against the 
peak flux density from the new Epoch HI observations. 
We plot peak flux density here instead of integrated flux 
density as variable sources are expected to be point-like, 
and the peak flux density is a better approximation of the 
actual flux density for unresolved sources. Each variable 
source is represented on this plot twice: once with its 
Epoch I flux density, and once with its Epoch II flux 
density. The solid line has a slope of log(x) = log(y) 
and should roughly bisect the sources if they are getting 
brighter or fainter at random. What we see, however, 
is a trend toward fainter Epoch HI flux densities. This 
is especially suspect, as the Epoch HI Stripe 82 observa- 
tions are the higher-resolution observations. To quantify 
this bias, we fit a line to the data, fixing the slope to 1.0, 
but allowing the intercept to vary. We found that the 
Epoch III observations have peak flux densities that are 
15% (±2%) fainter on average. The measured bias indi- 
cates that these sources may be resolved by the higher 
resolution Epoch HI data, necessitating the use of in- 
tegrated flux density for Epoch HI. Figure [2}d therefore 
shows the same sources, but with integrated flux density 
instead of peak flux density for Epoch HI. The bias we 
calculate (8% ± 2%) is now in the direction of slightly 
larger Epoch HI flux densities, and our estimate remains 
the same if we restrict the sample to those sources that 
are also point sources in Epoch HI. This offset is likely 
indicative of the systematic bias due to, for example, er- 
rors in absolute flux calibration. Due to the improvement 
in the correlation seen in Figure ^jp, we therefore find it 
necessary to use integrated flux density when comparing 
Epochs I/II with Epoch III. 

We then used the high-resolution Stripe 82 data to 
search for new variable sources. We matched the Epoch 
HI (Stripe 82) source catalog to both the Epoch I and 
Epoch II FIRST catalogs, again using a matching radius 
of 3". We required the sources to be point sources in 
the Stripe 82 data, where we defined a "point source" as 
anything with a ratio of peak-to-integrated flux density 
Siii,pk/Siii,int > 0.7. The exclusion of extended sources 
from a variability search is justified in that extended 
sources are not expected to be variable on such short 
timescales (less than 14/ [1 -I- z] years). Upon further in- 
vestigation, we found that this point-source criterion re- 
sults in the eli mination of some of t he previously-known 
variables (from Ide Vries et al.ll2004 ) from our candidate 
list due to the higher-resolution of the new data. It is 
therefore possible that some variable sources are missed 
by this cut. Indeed, we found that adjusting the value of 
this cutoff can have a large effect on the number of vari- 
able sources in our sample. However, for the purposes 
of this initial variability census, we argue that reliabil- 
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Fig. 2. — Peak (left) and integrated (right) flux density from the Epoch III St ripe 82 catalog versus peak flux density from FIRST (Epoch 
I or II) for the variable sources previously identified by Ide Vries et aTl II2004I ). Variable sources that appeared in both the Epoch I and 
Epoch II catalogs are represented twice. The typical error bars are much smaller than the plot symbols. The solid line indicates where the 
flux densities from the different epochs are equal. The dotted line is a fit to the data. The left panel shows that the Stripe 82 sources have 
peak flux densities that are 15% fainter on average than the corresponding peak flux densities from FIRST. In the right panel, a flt to the 
data is again overplotted as a dotted line, but it lies basically on top of the solid line and is not as noticeable. The bias indicated by the 
flt is 8%, this time in the direction of higher Epoch III flux density values. 
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Fig. 3. — Epoch III (Stripe 82) integrated flux density versus 
Epoch I or II FIRST peak flux density for all sources which satisfy 
our point source criterion and have a match between Epoch III and 
at least one of the earlier epochs. Sources in both FIRST epochs are 
plotted twice - once with each FIRST flux — and may be identified 
as variables with respect to Epoch III in one or both cases. Sources 
that exceed our preliminary variability threshold (Equation 1) are 
shown in pink , whil e sources previously-identified as variable by 
Ide Vries et al.l j 20041 ) are shown in cyan. 



ity is more important than completeness, and we 
required ah variable sources to meet this criterion. 

We defined a preliminary variable sample consisting of 
all point sources which satisfy the criterion: 



AS > 5 X (cTj^jj + a 



,1/2 



hi; 



(1) 



where AS is the flux density difference between the two 
epochs being compared, ctiji is the rms at the position 
of the source in cither Epoch I or II (depending on which 
is being matched), and crjn is the error on the integrated 
Stripe 82 flux density. Note that we also required a vari- 
ability amplitude of >5% to account for uncertainties in 
absolute flux calibration. For the FIRST epochs I and 
II, the rms is a good estimate of the error on the peak 
flux density, which we used as the best estimate of the 
true flux density. For Epoch III, on the other hand, we 
used the integrated flux density (as discussed earlier in 
this section), and since the uncertainty on the integrated 
flux density is greater than the rms listed in the catalog 
at the source position, we used the former as the flux 
density uncertainty. These errors were calculated by the 
automated source-finder H APPY during the c reation of 
the Stripe 82 catalog fsee iHodge et al.ll20lH) . For aU 
sources, the Poisson noise is negligable. 

We identified 1,436 distinct point sources (Figure [3]) 
with a match between the Epoch III catalog and at least 
one of the earlier epochs. Of these point source matches, 
258 satisfied the preliminary variability criterion defined 
in Equation 1. These variable source candidates (shown 
in pink) lie, by definition, on the outer envelopes of the 
distribution. One of the more obvious features of the 
plot is that the distribution of matches is asymmetric; 
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Fig. 4. — Top: A histogram showing the amphtude of the vari- 
abiUty (as a percent of the source strength) for the final variable 
sample. Due to the possible presence of systematic effects such 
as absolute flux calibration errors, we only consider sources which 
vary by more than 30%. Middle: The distribution on the sky of 
the final sample. The area with RA>350° is not covered by this 
study. Bottom: Distribution of the deconvolved major ax is values 
for th e previously-identified variable sources from de Vr ies et al.l 
1120041 ) (filled histogram) and the variable sources identified with 
the addition of Epoch III data (plain histogram). 



there are actually more sources that appear brighter in 
Epoch III. This is an effect of the Malmquist bias, and 
it was discussed in this context (Stripe 82 v ersus FIRST 
surve y) in the Stripe 82 VLA survey paper (jHodge et al.l 
1201 If ). For a source that is resolved in the Stripe 82 
data. Stripe 82 is not as sensitive as FIRST, which can 
cause the peak flux density of the source to fall below 
the detection threshold. It will appear in the Stripe 82 
catalog only if a noise fluctuation, or variability, causes 
the apparent peak flux to be higher. In that case, the 
source will also be brighter than the FIRST source. The 
bias seen in Figure [3] is therefore simply an artifact of the 
difference in resolution between the epochs. 

A second effect is visible at the bright end of the distri- 
bution. There it appears that bright sources tend to be 
fainter in the Stripe 82 epoch, even despite the use of in- 
tegrated Stripe 82 flux density. Due to the clear presence 
of resolution effects, therefore, we have conservatively de- 
fined as variable only those sources which have a higher 
flux density in the high-resolution (Epoch III) data. 

Finally, flux calibration plays an extremely important 
role in any study concerned with variability. If the abso- 
lute flux calibration differs between year to year or fleld 
to fleld, this could affect how many sources are flagged 
as variable. To estimate the effect of a flux calibration 
error, we tried increasing and decreasing the Stripe 82 in- 
tegrated and peak flux densities by small amounts. We 
found that if we altered the Stripe 82 fluxes by even 
5% (the value assumed above for the calibration uncer- 
tainty), the number of variable candidates changed dra- 
matically (by ~50%). If the absolute flux calibration is 
slightly different on certain fields, this might therefore 
cause a non-uniform distribution of the variable sources. 
We cannot attempt to correct for this effect by adjusting 
the fiux density ratios to a median value of unity since 
the the Malmquist effect and remaining resolution bias 
are uncertain. Because of this uncertainty, therefore, we 
imposed one final, additional constraint: we consider as 
variable only those sources with a flux density variation 
greater than 30%. Such large variations are less sensitive 
to smaller systematic effects like differences in absolute 
flux calibration. 

With these additional criteria, we defined a final vari- 
able sample of 89 sources. These 89 sources will be 
the variable sample we refer to for the remainder of the 
pa per. Of these s ources , 14 were previously identified 
by Ide Vries et al.l (J2004D . and 75 are newly-identified. 
Note that, because we only considered sources that are 
brighter in the high-resolution data, we do not find it 
necessary to disallow multiple Stripe 82 sources within 
the 3" radius, which could be due to a FIRST source 
being resolved into multiple components. If a FIRST 
source were to be resolved in the Stripe 82 data, the fiux 
density of the central component would appear fainter, 
not brighter. 

Some basic properties of the final variable sample are 
shown in Figured A histogram of percent variability dis- 
played by the variable sources is shown in Figure U] (top), 
where all sources vary by >30% as outlined above. The 
distribution on the sky of these sources is shown in Figure 
m (middle) . This plot confirms that the variable sources 
are not densely clustered in certain areas of the stripe, 
a result which would imply errors in photometric cali- 
bration rather than true variability. Figure |4] (bottom) 
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Fig. 5. — Epoch III (Stripe 82) integrated flux density versus 
Epoch I or II FIRST peak flux density for our flnal variable sam- 
ple. Sources from t he sample that were previously identified by 
Ide Vries et al.l I I2004I ) are shown in cyan. 

shows the distribution of deconvolved major axis val- 
ues. The open histogram shows the final variable sample, 
while the filled histogram shows the new ly-derived val- 
ues fo r the previously-identified variables (jde Vries et al.l 
l2004f) using the higher resolution Epoch III data. The 
distribution for our final variable sample extends up to 
major axis values of 2". Aside from the two outliers with 
Maj > 3", the previously-identified variables all have 
values of Maj < 1.7". We can therefore see that if we 
applie d our point source criterion to the Ide Vries et al.l 
(|2004l ) sample, we would be 96% complete, as 56 out of 
their 58 sources (covered here) have major axis values 
below 2". 

In summary, our final variable sample consists of 
sources which have a match (within 3") between the 
Stripe 82 data (Epoch III) and one or both FIRST 
epochs, and which meet the following four criteria: 

(1) Display fiux density differences greater than ba 
(Equation 1) 

(2) Point source in the high-resolution Epoch III data 
(i.e. Siii,pk/Siiiant < 0.7) 

(3) Brighter in the high-resolution Epoch III data (i.e. 
Siii > Slii) 

(4) Variability amplitude > 30% 

The application of these criteria results in 89 sources, 
and these source make up our final variable sample. 

3.2. Transient Phenomena 

Transient phenomena, although rare, can reveal inter- 
esting and new astrophysical sources when discovered. 
We searched for transient radio sources in our data in 
much the same way that we searched for variability. Due 
to the resolution bias, we only searched for sources that 
were undetected in the FIRST data, and not the other 
way around, since FIRST sources may simply be resolved 




1.4GHz Si„, (mJy) 

Fig. 6. — Variable fraction of radio sources as a function of flux 
density. The median flux density for each flux density bin is plot- 
ted. 

out of the Stripe 82 data. We considered only the Stripe 
82 sources that satisfied our point source criterion, and 
for each undetected source, we retrieved the FIRST rms 
in each epoch at that location from the respective cover- 
age maps. We required the source to satisfy our variabil- 
ity criterion, assuming an upper limit for the source of 
5xaFTRS T + 0.25mJyto account for the FIRST CLEAN 
bias (Bec ker et al.lll995l) . We also required source s to be 
absent from both FIRST epochs, as ide Vries et al.l (|2004D 
reported a systematic effect which caused a large number 
of sources to appear only in one FIRST catalog or the 
other at faint fiux densities. This left us with 14 initial 
transient source candidates. 

Upon close inspection of each of the candidates, none 
appeared to be real transient sources. The majority were 
sources that had been resolved in the Epoch III data such 
that the position of the FIRST component shifted by 
more than 3", or another component appeared 3" from 
the original component. This was not the case for two 
of the sources, which were clearly isolated point sources 
in the Stripe 82 images. However, we detected FIRST 
counterparts to both of these sources in the combined 
(Epoch I-I-II) FIRST catalog. Finally, one of the sources 
was a strong sidelobe {Sint = 4.3 mJy) of a bright source 
roughly 2' away. 

Therefore, we did not find any evidence for transient 
phenomena in the Stripe 82 data. We note that this re- 
sult is not particularly surprising, particularly since both 
the FIRST and Stripe 82 catalogs are based on coadded 
images using individual observations spread over days to 
months. Transients that appear on shorter timescales 
would therefore appear fainter in the coadded data, ren- 
dering them harder to detect. 

4. RESULTS 

Figure [5] is the main figure that characterizes the 
newly-defined variable sample. It shows the final vari- 
able sam ple, with those that were previously identified 
by de Vr ies et al.l ([MQ4) overplotted in cyan. We find 
that 14 of the 58 previously-identified variables satisfy 
our variability criterion. Note that this means they have 
increased in brightness, as we are limiting our variable 
sample to sources that are brighter in the high resolution 
epoch. Of the remaining de Vries et al. variable sources. 
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two sources did not survive our point source requirement 
and the remaining sources were either not found to be 
variable, or varied significantly in flux density but were 
fainter in Stripe 82, and thus ambiguous. Assuming that 
the number of sources that get brighter is roughly equal 
to the number of sources that fade with time, we esti- 
mate that the true number of variable sources (based on 
our other criteria) is double the number we find here, or 
approximately 28 sources of the 58 previously identified. 

Table [T] lists the 89 variable sources discovered by 
matching the two FIRST epochs to the Epoch III data. 
Variable sour c es th at were previously discovered by 
Ide Vries et al.l ()2004[ ) are indicated with an asterisk be- 
fore the RA. Values for the peak and integrated flux den- 
sity and rms for each of the three epochs are listed, as 
well as the deconvolved major axis value measured from 
the Epoch III data, the variability amplitude /var (Sec- 
tion |42]), and the morphological classification in FIRST 
and Stripe 82 (Section 14. 5p . Values of 0.0 indicate that 
the source was not detected in that epoch. Note that 
strong sources significantly affect the local noise, causing 
the rms listed to increase dramatically from the median 
survey value of 52 /iJy. Assuming that an equal num- 
ber of sources increase and decrease in flux density, we 
would expect approximately 178 variable sources total in 
the area, or ~3 variables deg~^. 

Note that we used a 5a var iability threshold instead 
of the 4cr threshold used by de Vries et all (l2004f) . (They 
also used a slightly different equation - see Ide Vries et all 
120041 for details.) However, Epoch III has a lower typ- 
ical rms, and is therefore sensitive to lower levels of 
variability. This explains why a small number of the 
newly detected variable sources were actually slightly 
more variable between Epoch I and II than they were 
be tween Epochs I/II an d III, and yet were not identified 
bv lde Vries et al.l ()2004[ ) as variable sources. 

T able [2] lists all of the 58 sources that were part of 
the Ide Vries et al.l ()2004[ ) sample and occur in the area 
covered by the new Epoch III data. The columns are 
the same as Table [TJ Sources that meet our variability 
criteria for Epoch III (including a brighter flux density 
in the new epoch) are again indicated with an asterisk 
before the RA. As the sources in this sample tend to be 
much brighter, in general, than the sources selected with 
the Epoch III data, the local rms values are also higher. 

4.1. Dependence on Source Strength 

To investigate whether the variable fraction varies with 
flux density, we split the variable sample into four flux 
density bins. The median flux values of the bins ranged 
from 4 mjy on the faint end to 250 mJy on the bright end. 
The variable fraction for each bin is plotted against the 
median flux density in Figure[6l We find that the variable 
fraction remains essentially constant as a function of flux 
density, with a possible increase in the fraction of variable 
sources at fainter flux densities (though we caution that 
a slightly different binning can remove this subtle effect, 
and it may therefore not be significant). This is despite 
the fact that the data are inherently more sensitive to 
variability in very bright sources, while for faint sources 
it takes a large change in amplitude to achieve the 5cr 
threshold. Our requirement of 30% variability makes the 
sample more consistent, but the variable fraction is very 
sensitive to the exact cut applied. Therefore, it appears 
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Fig. 7. — Histogram of fractional variability of the variable 
sources, defined as the maximum integrated flux density over the 
minimum integrated flux density. Note that since we require 
sources to vary by at least 30%, nothing has a fractional variability 
of less than 1.3. 

that faint radio sources are just as variable as bright radio 
sources, and possibly more so. 

4.2. Fractional Variability 

A paper by lBecker et al.l (|2010( ) on variable sources in 
the Galactic plane found that Galactic sources were more 
variable than extragalactic sources. Their Galactic sam- 
ple came from three epochs of VLA data on the Galactic 
Plane taken roughly 15 years apart, and they co mpared 
this against the extragalactic sample of de Vries et al.l 
(JOOjl;) . To see if this still holds true for our more robust 
sample of extragalactic variables, we calculated the frac- 
tional variability for the sources in our variable sample, 
where fractional variability is defined as: 



Jva 



(2) 



where S = Spk for the FIRST epochs I/II and S = Si„t 
for Epoch III. 

Figure [7] shows the distribution of values for the frac- 
tional variability. Note that since we required a vari- 
ability amplitude of > 30% to be in the variable sam- 
ple, no sources have a fractional variability less than 
1.3. Our values are c ompared to the extragalactic sam- 
ple of Ide Vries et al.l (|2004') and the galactic sample of 
IBecker et al.l (I2010D in Table El T he bins are as they 
were defined in IBecker et al.l (|2010f) . The lowest bin is 
empty for our sample and the second-lowest bin only in- 
cludes sources with fractional variabilities down to 1.3 
due to our variability amplitude requirement. Never- 
theless, there is a clear difference in the distribution of 
fractional variability between our sample and the Galac- 
tic sample. Ignoring the lowest bin, almost half (46%) 
of our extragalactic sources have a fractional variability 
fvar < 1.5, while only 13% of Galactic sources fall in the 
same range. Meanwhile, only 1% of extragalactic sources 
have a fractional variability f^ar > 3 (corresponding to a 
variability amplitude of 200%), compared to a full 44% of 
Galactic sources. The difference between the samples is 
significant at the >5.5a- level based on a x^ contingency 
test. Therefore, we confirm that Galactic sources appear 
to be more highly variable than extragalactic sources. 
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Fig. 8. — Mean fractional variability as a function of flux density. 
We took the average flux density for each individual variable source, 
and the plotted points show the average of those individual values 
within a bin. 

We show the flux dependence of the fractional variabil- 
ity in Figure [H where we plot mean fractional variability 
as a function of flux density bin. At low flux densities, 
where most of the variable sources lie, there appears to be 
a steep rise in fractional variability. Since our variable 
sample has roughly the same flux density distribution 
as the Galactic sample, this trend cannot explain the 
very different fractional variability distributions of the 
two samples. Within our sample, it is likely a selection 
effect, since (as discussed in Section |4TT|) fainter sources 
must exhibit a larger (fractional) change in amplitude in 
order to exceed the 5cr threshold in the flrst place. 

4.3. Characteristic Timescale 

We can attempt to infer something about the char- 
acteristic timescale of AGN variability by looking at the 
number of variables discovered between the three epochs. 
Splitting our results up by epoch, we flnd 55 variable 
sources between Epochs II and III (a 7 year baseline) 
and 74 variable sources between Epochs I and III (a 14 
year baseline). These numbers do not add up to the 89 
sources in our final sample because some sources regis- 
tered as variable in both searches. 

However, we cannot directl y compare the s e nurn bers 
to the 58 variables found bv Ide Vries et al] (|2004[ 1 be- 
tween Epochs I and II since we used different variability 
criteria and the new data reach a lower rms sensitivity. 
For a fairer comparison, we redid the search for variables 
between the FIRST epochs (I and II) using our Equa- 
tion 1 and requiring >22% variability. (Note that the 
requirement for comparisons with Epoch HI was 30%, 
but the small, 8% bias in the direction of higher Epoch 
HI fluxes (Figure [2]) means that the correct value to use 
here is reaUy 30% - 8% = 22%.) To address the differ- 
ent rms sensitivities reached, we also redid our variability 
search between Epochs I/II and Epoch HI, substituting 
the FIRST rms of 0.15 mJy for the Epoch HI rms and 
scaling it up by VNbcams since we are using the integrated 
Stripe 82 flux. Finally, we multiplied these results by a 
factor of two, assuming that the same number of sources 
fade as get brighter. 

With these changes, we flnd 64 variable sources be- 
tween Epochs I and II, 80 variables between Epochs II 



and HI, and 104 variables between Epochs I and HI. The 
timescale between Epochs I-II and II-III is seven years 
in both cases, so we would expect to find roughly the 
same number of variables each time. The numbers we get 
are consistent within 25%, with the remaining difference 
likely explained by the intricacies of selecting a consistent 
sample across datasets with different angular resolutions 
and sensitivity limits. The number of variables discov- 
ered between Epochs II and III (80) and Epochs I and 
HI (104) are more reliably compared, as they both in- 
volve a FIRST epoch matched to the Stripe 82 epoch, 
and the flux calibration of the two FIRST epochs has 
been verifled against one another (requiring only a small 
correction factor - see Section 2). Of the 104 sources 
discovered using the 14-year time span, 62 were already 
detected in the 7-year span, and 42 are newly detected. 
It thus appears that lengthening the observing timescale 
by 7 years only produced <50% more variables. While 
we do not have enough information to be quantitative, it 
thus appears that the average characteristic timescale for 
variability for these sources may be less than 14 years. 

4.4. Optical Counterparts 

We now turn to the optical properties of the sam- 
ple. To investigate the nature of the variable sources, 
we matched our sample to three different optical cat- 
alogs: the general catalog of all sources from the 
coadded imaging data in Stripe 82, a catalog of 
spectroscopically-confirmed quasars from SDSS, and a 
catalog of photometrically-selected quasars from SDSS. 
In order to interpret the results, we created a control 
sample of non-variable sources with which to compare. 
As with the variable source sample, we required the ra- 
tio of peak to integrated flux density Speak / •S'int > 0.7. 
Whereas for the variable sample we required the change 
in amplitude to be > 5(7, for the non-variable sample we 
required the change in amplitude to be < 0.5i7. Since 
variability is really a spectrum, this restriction ensures 
that we are comparing to the least variable objects pos- 
sible. 

We further required that the control sample have the 
same flux density distribution as the variable sample. To 
do this, we binned the integrated flux density of the vari- 
able sample into logarithmic bins of 0.3 dex distributed 
from 1 mJy to 1000 mJy, and we required the control 
sample to follow the same distribution. We used a ran- 
dom number generator to select the non-variable sources 
from each bin to go into the control sample. 

If the variable sources are brighter than the control 
sample sources in the optical, this could bias the re- 
sults, causing a higher overall match rate. In addition, 
brighter sources are easier to classify into galaxy/stellar 
sources, with fainter sources more likely to be classifled 
as galaxies. To see if the magnitude distributions for the 
variable/non-variable samples are similar, we retrieved 
SDSS r-band magnitudes for the sources with optical 
matches in the variable and non-variable samples. The 
cumulative distribution function of these magnitudes is 
shown in Figure [9] for both the variable sources (solid 
line) and non-variable control sample (dashed line). We 
used the Kolmogorov-Smirnov test to determine if the 
two magnitude distributions were likely to be drawn 
from the same underlying distribution. The resulting KS 
statistic is 0.12 with a probability of 0.85 (where small 
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Fig. 9. — Cumulative distribution function of SDSS r-band mag- 
nitude for the variable sources (solid line) and non-variable control 
sample (dashed line). 

values of the probability indicate that the cumulative dis- 
tribution functions are significantly different). We con- 
clude that the difference between the magnitude distribu- 
tions is most likely not significant and should not signif- 
icantly influence star/galaxy classification or the overall 
match rate. 

We first matched the variable and control samples to 
the optical Stripe 82 catalog, created from combining the 
multiple epochs of SDSS imaging data of Stripe 82 into a 
single catalog that extends tw o magnitudes deep er than 
the main SDSS imaging area (jAnnis et al.ll201ll) . After 
trimming this deep catalog to the area covered by the 
variable source sample, there are still over four million 
sources in the resulting optical catalog. By matching to 
these sources, we find that the overall match rate for the 
variable source sample is 60% ± 5%, while the match rate 
for the control sample is 44% ± 5%. Since it has been 
shown that the astrometry of the Stripe 82 radio catalog 
is well-tied to the SDSS (Hodge et al. 20lT|), we used a 
matching radius of 1" to minimize the number of false 
(coincidence) matches. To determine the rate of coinci- 
dence matches, we generated a fake catalog by shifting 
the radio source positions by 1' in four different directions 
and repeating the matching, averaging the results. The 
false match rate we determined is 2%, meaning that the 
vast majority of the matches reported above are likely to 
be real. 

To see if the results depend on flux density, we plot 
SDSS match fraction as a function of flux density in Fig- 
ure 1101 The variable sources show an increase in match 
fraction for brighter sources, while the non- variable 
sources show a matched fraction that is independent of 
radio flux density. While the matched fraction for the 
two samples appears consistent at low flux densities, the 
variable sources clearly have a higher matched fraction 
than the non-variables at high flux densities. 

The SDSS catalog includes a classification for each 
source based on star-galaxy separation, which differen- 
tiates between resolved and unresolved sources using the 
difference between the PSF magnitudes and model mag- 
nitudes (either de Vaucouleurs or exponential). For the 
radio sources that were optically identified, we utilized 
the SDSS star-galaxy separation to learn more about 
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Fig. 10. — Fraction of sources that match to the SDSS Stripe 82 
optical catalog as a function of integrated flux density. Variable 
sources are shown by the solid line, and the non— variable control 
sample are shown by the dashed line. The data were split into two 
bins, and the data points indicate the centers of those bins. 
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Fig. 11. — Fraction of optically-detected radio sources that are 
classified as stellar based on the SDSS star/galaxy separation. 
Variable sources are shown by the solid line, and the non— variable 
control sample are shown by the dashed line. The data points again 
indicate the centers of the bins. 

the optical counterparts, taking sources classified as 'stel- 
lar' to be quasar candidates. The results are as follows: 
Of the variable sources with optical counterparts, 72% 
are classified as galaxies and 28% are classified as stellar 
based on the star-galaxy separation, with a formal error 
of ±6%. For the control sample, we find that the identi- 
fied sources are composed of 92% galaxies and 8% stellar 
sources (±4%). These results are displayed along with 
the SDSS match fractions in Table HI 

These results demonstrate that variable sources are 
more likely to be quasars than non-variable sources. 
They also show that, despite this fact, a significant ma- 
jority of the variable sources appear to be extended in the 
optical. This could indicate that these sources are AGN 
in galaxies as opposed to quasars. We should caution 
that for fainter sources, the SDSS star-galaxy separation 
is not as reliable as for bright sources. In addition, the 
distinction may not be so clear-cut for these deep opti- 
cal data. Many low-z quasars that were point sources in 
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the single epoch data are now extended sources, because 
the stacking brings out the galaxy. However, we can still 
learn something about the underlying source population 
by looking at the flux dependence with respect to the 
control sample. To see if the stellar fraction depends 
on flux density, we plot stellar fraction as a function of 
flux density in Figure[TT] We split the data into two bins, 
and the plotted points represent the centers of those bins. 
Although the stellar fraction increases with flux density, 
it is below 40% for both bins. We therefore find that 
optically-extended sources make up the majority of vari- 
able sources, regardless of flux density, but a significant 
minority of sources are still quasars that outshine the 
host galaxy. We also see that the stellar fraction of the 
variable sources is larger than that of the non-variable 
sources for both high and low flux densities, although 
the error bars are large enough that this result is not 
statistically significant. 

To explore the quasar fraction further, we also 
matched the variable and control samples to a sample 
of spectroscopically-identified quasars from the SDSS. 
We used the fifth editio n of the SDSS quasar cata- 
log ([Schneider et al.ir2010[ ). which is based on the Sev- 
enth Data Release (DR7). The catalog contains 2083 
spectroscopically-confirmed quasars in the area covered 
by the variable source sample. We find that 9 out of 
89 variable radio sources match to a spectroscopically- 
confirmed quasar within 1", or 10% ± 3% (Table [5]). For 
the non-variable control sample, the percentage is 3% 
± 2%. Although limited by small-number statistics, we 
find that the variable sources in our sample are, again, 
more likely to be quasars than the radio sources in the 
control sample. 

As a final comparison, we matched the variable and 
non-variable control samples to a larger catalog of 
photometrically-selected quasars. The SDSS sample of 
p hotometrically-selected quasars comes from the work 
of [Richards et all (|2009D . The catalog is expected to be 
~80% efficient and contain 850,000 bona fide quasars in 
the area covered by SDSS DR6. There are approximately 
7,000 photometrically-selected quasar candidates in the 
area covered by our observations. Matching against this 
catalog, we find that 13% ± 4% of the variable sources 
correlate with photometrically-selected quasars, versus 
2% ±2% of the control sample sources. Here, the quasar 
fraction of the variables is (again) higher than that of 
the control sample (though we caution again that we 
are dealing with small number statistics). The control 
sample actually matches to the photometric catalog at a 
lower rate than it matches to the spectroscopic catalog 
even though the photometric catalog contains many of 
the spectroscopically-confirmed quasars. This is because 
the spectroscopic catalog contains some quasars that are 
not in the photometric catalog, presumably because they 
exhibit abnormal colors for quasars and were targeted 
for spectra for some other reason. The control sample 
matched to several of these quasars, while not matching 
to any purely photometrically-selected quasars, resulting 
in a lower matched fraction than when it was matched 
against the spectroscopic catalog. 

Perhaps a more interesting observation is that, despite 
having 3.5 times the number of quasar candidates as 
the spectroscopic catalog, the photometric quasar cata- 
log produces only ~50% more matches with the variable 



sampl e. This is the exact same effect seen bv lHodge et al.l 
(|2011f ) when matching the entire Stripe 82 catalog to 
both the spectroscopic and photometric quasar catalogs. 
There, it was shown that the two samples are actually 
entirely consistent to roughly r = 20, meaning that the 
discrepancy lies entirely with faint sources. The likely 
explanation has to do again with the SDSS targeting 
pipehne. The SDSS targets sources for spectra for a 
number of reasons, and the selection is quite efficient. 
This leaves the photometric catalog to fill in mainly those 
objects without radio detections, which are presumably 
more numerous. 

4.5. Morphology 

The higher angular resolution of the Epoch HI Stripe 
82 data means that we can investigate in more de- 
tail the morphology of both the newly-identified and 
the previously-known variable sources. In order to 
characterize t h e mor phology of the sources, we follow 
Ide Vries et al.l (|2004f ) and define the following six mor- 
phological classes: (PS) - Isolated point source; (CJ) 
- "core-jet" morphology, either as two separate com- 
ponents or as a single component with a clear elonga- 
tion; (CL) - "core-lobe" morphology, where the core 
is surrounded by two distinct lobe components, which 
are not variable; (CH) - "core-halo" morphology, con- 
sisting of a point source core surrounded by a diffuse 
radio halo; (CX) - complex source, consisting of mul- 
tiple components; and lastly, (HS) - hot-spot variabil- 
ity. These morphological classes are indicated for the 
new and previously-identified variable sources in Tables 
[Hand [3 

Of the 58 previously-known variables in the area cov- 
ered by the new Stripe 82 radio data, 45 were classified 
(based on the FIRST imaging) as having PS morphol- 
ogy, and 13 fell into other classes (8 CJ, 4 CL, and 1 CX; 
Ide Vries et al.|[2b04f ). With the higher-resolution Stripe 
82 data, we resolve four of the C J sources enough to now 
see a CL morphology, as well as one of the CL sources 
into a CX morphology. For three of the sources previ- 
ously classified as PS, we resolve the radio source into 
CJ morphology. Only one source moved into a simpler 
morphological class, and that is the source previously 
classified as CX, which now appears to be a PS with an 
unrelated neighbor. The breakdown we find is therefore 
42 PS (70.7%), 7 CJ (17.2%), 7 CL (10.3%), 1 CX (1.7%), 
and no CH or HS. One source was not detected in the 
Stripe 82 catalog and therefore was too faint to classify. 
All of the variability is consistent with being due to the 
radio core/AGN component. 

We next look at the morphology of the 89 variable 
sources identified in this work by the addition of a third 
epoch of r adio data. Note that 14 were identified pre- 
viously bv Ide Vries et al.l ([20041 ). but 75 are newly dis- 
covered. We find that four sources fall into the CJ class 
(~5%), and five quahfy as CL (~5%). No sources show 
any of the more unusual morphologies (CH, CX, or HS). 
Therefore, a full 90% of our variable sample are simple 
point sources. 

Figure [T2| shows contour plots for each of the 
previously-known and newly-identified variable sources 
that have non-PS morphology. The images are taken 
from the high resolution Stripe 82 radio data. The 
position of the variable component is indicated with a 
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Fig. 12. — Contour plots of the variable sources with complex morphology, using data from the Stripe 82 survey. The cross indicates the 
position of the variable component . Contours are 60 /ijy X ±(3,5,7,9,12,15,20,30,50,100,300,1000), where black contours are positive and 
grey contours are negative. 
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Fig. 1121 (continued) 
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Fig. 13. — Greyscale image (zoomed-in) of complex source 
J222729+000522. 

cross. Some of the images (e.g., J012213-001801 and 
J220755-00215) exhibit strong negative (and positive) 
contours characteristic of incompletely-cleaned sidelobes 
and missing short spacings. These features are clearly 
image artifacts and were accounted for in the morpholo- 
goical classification. Further notes on individual variable 
sources, including the complex source J222729+000522 
(Figure llSp can be found in the Appendix. Any ref- 
erences to optical matchin g refe r to the deep coadded 
SDSS catalog ()Annis et al.ll201lD . 

5. DISCUSSION 

Our study of lon g-term AGN variability extends the 
work of lde Vries et al. (2004) by introducing a new epoch 
of high-quality, high-resolution data. One thing that 
the addition of a third epoch allows us to do is to 
constrain the characteristic timescale of variability for 
these sources, providing insight into the mechanism for 
variability as well as helping to identify the type of 
source. For example, there is some evidence that QSO 
have longer characteristic timescales than BL Lac ob- 
jects (jHughes et al.llT992l ). Because we are measuring 
timescales in our reference frame instead of transform- 
ing to the source frame, we are introducing an additional 
source of dispersion into the data. However, no known 
correlation exists between redshift and timescale beyond 
1 + z, and intrinsic beaming is likely to be a more impor- 
tant effect. 

We find that extending the observing span from 7 to 
14 years produces <50% more variables. From this, we 
can infer an upper limit on the characteristic timescale 
for AGN variability of 14 years. We compare this re- 
sul t to the long-term monitoring campaign carried out 
bv lAUer et al.l (I1992D on the Univ ersity of Michigan 26m 
paraboloid. iHughes et al.l ()1992f ) performed a structure 
function analysis on the QSO and BL Lac objects in the 
sample, finding mean characteristic timescales of 1.95 
and 2.35 years, for BL Lacs and QSOs, respectively, 
though with very broad distributions. They found that 
the majority of sources with timescales greater than 10 
years were QSOs . In a longer-term analysis of the same 
data, lAUer et all ()2003D reported that the variability of 
the galaxies in the sample consisted of slower variations. 



often only identifiable over a few decades. They at- 
tributed this to relatively long characteristic timescales 
in combination with low luminosity of the AGN relative 
to that of the extended structure, requiring unusually 
bright events for detection. Both papers were based on 
a very bright (55ghz > 1.3 Jy) sample of objects, how- 
ever. Our study is the first study to try to determine the 
characteristic sample of such a faint (and blind) sample. 
Additional epochs are required in order to constrain the 
characteristic timescale further. 

The fraction of sources that are variable at the milljan- 
sky level is another unknown. At high flux density levels, 
variable sources typically correspond to flat-spectrum 
AGN such as blazars and flat-spectrum QSOs. How- 
ever, there is a general consensus that the nature of 
the radio source population changes with decreasing flux 
densit y. The source cou nts flatten below a few millijan- 
skies (|Windhorst et all [19851 . which is generally taken 
as evidence that a new type of source dominates be- 
low that level. The exact composition of the popula- 
tion and ho w it depends on flux density is still hotly de- 
bated (e.g IStrazzullo et al.l 120101: IPadovani et"al] 120091: 
llbar et aLl 120091 ). What is clear is that around a flux 
density of ~1 mJy, flat-spectrum AGN make up only a 
small percentage of the population, with some mixture of 
star-forming galaxies and steep-spectrum sources mak- 
ing u p the remainder fe.g. lHopkins et aLll2000t iPadovanil 
120111 ). One might therefore naively expect the variable 
fraction of sources to decrease with flux density as well. 
Our results, on the other hand, imply that the fraction 
of sources that are variable stays constant down to low 
fl ux densities, and po s sibly even increases slightly. 

iThvagaraian et all ()2011l ) have published the largest 
sample of millijansky variable sources to date, but they 
do not investigate how the variable fraction depends on 
flux density. The only other blind variability survey of 
sufflcient size a.nd d epth to compare to is the work of 
Ide Vries et all (|2004[ ) on the two FIRST epochs. How- 
ever, they specifically adjusted their selection criteria to 
force the variable fraction to be constant, so a compar- 
ison is out of the question. This is therefore the first 
large-scale study (to our knowledge) investigating the 
dependence of radio variability on flux density down to 
the millijansky level. We conclude that the variable frac- 
tion does not decrease at faint flux densities despite the 
changing composition of radio sources. This may im- 
ply that the variable fraction will remain constant (or 
even increase) down into the sub-millijansky regime for 
next generation surveys l ike that planned by ASKAP and 
MeerKAT (although c.f. lCarilh et~al|[2003f) . 

In addition to variability, the existence of multiple 
epochs of data allows us to look for the more rare, tran- 
sient events as well. Because of the resolution bias, we 
were only able to search in the Stripe 82 epoch, but we 
found no evidence for any such events. Ide Vries et al.l 
(|2004l ) also reported finding no evidence for transients 
down to 2 mJy in their search over the larger area cov- 
ered by bot h FIRST epoc hs. 

Recently, iBecker et al.l (j2010( ) showed that Galactic 
sources are much more highly variable than extragalactic 
sources. We have confirmed this result, finding that only 
1% of extragalactic sources have a fractional variability 
>3, while 44% of Galactic sources vary by this much. 
The distribution of variability amplitude seen in our sam- 
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pi e is roughly co nsistent with the extragalactic sample 
of Ide Vries et al.l ((200 4.1 used in the original comparison, 
although we see slightly more variability of intermediate 
amplitu de. We cons ider our results to be more reliable 
than the Ide Vries et alJ p004l ) results, as our images are 
much higher quality than the FIRST images used exclu- 
sively there. 

We saw in Section 14.41 that the overall optical identi- 
fication rate of the variable sources is higher than that 
of the non-variable control sample, and that the effect 
increases at higher flux densities. Our match rates for 
the h igher flux dens i ty bin agree with the overal match 
rates Ide Vries et al.l (|2004[ ) found for their brighter (on 
average) variable sample - they reported that 82% of the 
variables were identified, versus only 38% of the control 
sample. Our results imply that our variable sample con- 
tains a higher fraction of quasars than our non-variable 
sample, particularly at higher flux densities. As quasars 
are unresolved, an imaging survey will detect them more 
easily than extended sources of the same (integrated) 
magnitude. As the unidentified radio sources are most 
likely obscured AGN, our results also imply that variable 
sources are less likely to be obscured AGN than non- 
variable sources, though the effect decreases at lower flux 
densities. 

Of the optically-identifiable variable sources, we find 
that a majority match to extended sources (galaxies), 
as opposed to point sources (qua sars). In cont r ast, a 
full half of the variables studied bv lde Vries et al.l ()2004l ) 
were optical point sources. Our sample is fainter than 
the de Vries et al. sample in the radio, on average, 
which could explain some of the discrepany. Not only 
are fainter sources more likely to be galaxies, but the 
SDSS star-galaxy separation is less reliable at low flux 
densities. However, when we calculate the stellar frac- 
tion as a function of flux d ensity, it never appro aches the 
higher percentages seen in Ide Vries et al.l ()2004l ) even for 
high values of flux density. This may be largely because 
the coadded optical data we compare to go much deeper, 
allowing us to identify more of the counterparts (not just 
the bright quasars), but also allowing the host galaxies 
of many low-z quasars to be detected, and blurring the 
line between traditionally-defined "quasars" and generic 
AGN in galaxies. 

We also compare these results to the optical point 
source fraction of a control sample, defined to be non- 
variable in the radio. We find that the variable sources 
are more likely to match to point sources (quasars) than 
control sample sources. This is further confirmed by 
matching both samples to actual spectroscopic and pho- 
tometric quasar catalogs (though with the added caveat 
that the results suffer from small number statisti cs). 
This result is consistent with the results of de Vri es et aD 
(J20041 ■ who also used the SDSS star-galaxy separation 
to classify their identifiable fraction and reported a sig- 
nificantly higher stellar fraction in the variables than in 
the control sample. 

Upon examining the morphology of the previously- 
identified variable sources using the new, higher- 
resolution data, we find that 26% have structure that 
is more complex than a simple point source. We identify 
new jet/lobe components in a number of the variables, 
resulting in an increased fraction of complex sources over 
the 21.8% reported by ,dc Vries et al, (,2004 ). This con- 



firms their flnding of a strong correlation between radio 
variability and jet activity for this sample of radio vari- 
ables. Furthermore, as we resolve four of the sources 
previously classifled as CJ enough to see radio lobes, we 
believe that deeper data of the same angular resolution 
would flnd than an even higher fraction of variables ac- 
tually reside in the cores of FR II radio galaxies. 

We also examine the morphology of our sample of 
89 variable sources. According to the unified model 
of AGN, extended emission would arise from the radio 
lobes, which are not expected to be variable. All of the 
variability we see is consistent with originating from the 
core, which is consistent with the model. A full 90% of 
the variables are, in fact, simple point sources, as op- 
posed to being cores of FRIIs, or showing jet-like fea- 
tures. This is significantly different from the results for 
the previously-identifled variables, and may be due to a 
couple of reasons. For one thing, our point source crite- 
rion results in a stricter cut due to our higher resolution, 
and we are therefore more biased against sources which 
are extended due to the presence of jets. The higher an- 
gular resolution data may also resolve out some of the 
extended structure. Finally, our sample is fainter, on 
average, than the de Vries et al. sample. Although the 
Epoch III data are three times more sensitive, the fac- 
tor of three increase in angular resolution means that 
diffuse components are spread out over nine times the 
area. This would make diffuse components of the same 
fractional strength much harder to detect. 

6. CONCLUSIONS 

We carried out a blind survey for extragalactic ra- 
dio variability at the millijansky level by comparing two 
epochs of data from the FIRST survey with a new 1.4 
GHz survey of SDSS Stripe 82. The three epochs, which 
we refer to as I, II and HI, are spaced seven years apart 
and cover an overlapping area of 60 deg^. Our main re- 
sults are the following: 

• We identifled 89 variable sources with a variability 
amplitude above 30%. Of these, 14 s o urces were 
previously reported by Ide Vries et al.l ()2004( ) and 
75 are newly-identified. Due to the resolution bias, 
our sample is restricted to those sources that in- 
creased in flux density in Epoch HI. Assuming that 
an equal number of sources increase and decrease 
in fiux density, we would expect approximately 178 
variable sources in total, or ~3 variables deg~^. 

• We found no evidence for transient phenomena in 
the new Stripe 82 data. 

• The fraction of radio sources that are variable ap- 
pears to remain constant down to the millijansky 
level, possibly even increasing slightly at low flux 
densities. This is despite the fact that the fraction 
of flat -spectrum (i.e. canonically variable) sources 
decreases with decreasing fiux density. This may 
imply that next generation radio surveys with tele- 
scopes like ASKAP and MeerKAT will see a con- 
stant or even increasing variable fraction down into 
the submillijansky regime. 

• By comparing the distribution of fractional vari- 
ability against that of a Galactic sample, we con- 
firmed that Galactic sources tend to be more highly 
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variable than extragalactic sources. While 44% of 
Galactic sources vary by more than 200%, only 1% 
of extragalactic sources have variability amplitudes 
this large. 

• When we compared the number of variables dis- 
covered over different lengths of time, we found 
that increasing the observing interval from 7 to 14 
years only resulted in the detection of <50% more 
sources. From this, we inferred an upper limit to 
the mean characteristic time of AGN variability of 
14 years. 

• We found that the overall rate of SDSS optical iden- 
tification is higher for the variable sources than 
for a flux-matched non-variable control sample, 
though the effect depends on flux density. This 
implies that the variable sample contains a higher 
fraction of quasars than the non-variable sample, 
a result which is confirmed by both the SDSS star- 
galaxy classification and comparisons to spectro- 
scopic and photometric quasar catalogs. The ma- 
jority of variable sources at all flux densities appear 
to be extended galaxies in the optical, implying ei- 
ther that quasars are not the dominant contributor 
to the variability of this sample, or that the deep 
optical data allow us to detect the host galaxies 
of many low-z quasars, blurring the line between 
traditionally-defined "quasars" and generic AGN. 

• We used the new epoch of high-resolution data to 
characterize the radio morphology of the variable 
sample, concluding that over 90% of our variable 



sample are stand-alone point sources. This sample 
does not appear to exhibit the strong correlation 
with radio jets seen in the only other blind survey 
of millijansky AGN variability. The resolution bias 
is likely partly to blame. 

While this work serves as a starting point for the study 
of faint AGN variability, clearly more data are needed 
if we wish to better understand the variable population 
and characterize the variability. Deeper, high resolution 
data are needed to identify the population, and addi- 
tional epochs are essential to constrain the characteristic 
timescale and mechanism for variability. Nevertheless, 
this survey gives us a small taste of what is to come 
with variability campaigns on upcoming telescopes like 
ASKAP and MeerKAT. 
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APPENDIX 



NOTES ON INDIVIDUAL VARIABLE SOURCES 



J005205+003538 (CL): Previously known to be variable. Diffuse lobes are now visible on either side of the core. 
Although neither component appears in the Stripe 82 catalog, we believe this fits better in the CL category than the 
CJ category, as two diffuse lobes are clearly visible. The central (variable) component corresponds to an optical QSO 
at z = 0.399 with a magnitude of r = 16.3. 

J005717-002433 (CJ): Previously known to be variable. The contour plot shows that the FIRST source has been 
resolved into three components, although only the central source shows up in the catalog. There is a core, a bright 
point source connected to the NE, and a separate extended component to the SW. Because the relationship of the 
components is unclear, we classify this as a core-jet system. The central (variable) source matches to a QSO at 
z = 2.712. 

J010155+004129 (CJ): A newly-identified variable source. The higher-resolution Stripe 82 data show a diffuse, 
jet-like stream connecting the central source to a hotspot in the north. This stream consists of three components in 
the Stripe 82 catalog: a component 7" away with a peak flux density of 0.4 mJy beam~^ and an integrated flux density 
of 2.3 mJy; an elongated component 10" away with a peak flux density of 0.4 mJy beam~^ and an integrated flux 
density of 6.3 mJy; and a third component 12" away with a peak flux density of 0.6 mJy beam"^ and an integrated 
flux density of 2.6 mJy. The central (variable) radio component corresponds to an optical QSO at z = 0.649. 

J012213-001801 (CL): Previously known to be variable. The FIRST image showed a core and two lobes, one (the 
northeastern lobe) much fainter than the other (peak flux density of 1.7 mJy beam~^). In the Stripe 82 data, the 
northeastern lobe is (mostly) resolved out and fainter than the nearby sidelobes. It does not make it into the Stripe 
82 catalog. The central (variable) component has a peak flux density of 331.3 mJy beam~^ and an integrated flux 
density of 405.0 mJy. It coincides with a photometrically-selected z = 1.325 UV-excess (UVX) quasar candidate 
([Richards et al.|[2004[ ). The southwest component has a peak flux density of 38.2 mJy beam~^ and an integrated flux 
density of 116.2 mJy and does not have an optical counterpart. 

J012956+002338 (CL): A newly-identified variable source. This radio-double has an FRI morphology, with the 
lobes terminating in hot spots and a faint jet-like feature visible near the north-eastern lobe. The north-eastern lobe 
is visible in the image, but has been resolved out somewhat and does not appear in the catalog. The south-western 
lobe has been resolved into three distinct components: a component 19" away with a peak flux density of 1.0 mJy 
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beain~^ and an integrated flux density of 1.9 mJy; a component 22" away with a peak flux density of 0.8 mJy beam~^ 
and an integrated flux density of 1.1 mJy; and a third component 23" away with a peak flux density of 1.2 mJy beam~^ 
and an integrated flux density of 2.1 mJy. The variable component (core) corresponds to a QSO at z = 1.079. 

J013100+002747 (CL): A newly-identified variable source. The FIRST data show only the central point-source, 
but the Stripe 82 data show evidence of a double-lobed structure, although the southeastern component is only 
detected at la significance. The central (variable) component matches (at 0.07") to a galaxy of magnitude r = 21.7 in 
the deep optical catalog. The northwestern component has a peak flux density of 0.8 mJy beam~^ and an integrated 
flux density of 1.4 mJy. This component matches (within 0.04") to a stellar source of r = 24.2 in the deep optical 
catalog, suggesting that it may be unrelated. In this case, the morphological classification would be PS. 

J013814+001444 (CJ): Previously known to be variable. The FIRST source has been resolved into two components. 
The central (variable) component has a peak flux density of 18.9 mJy beam~^ and an integrated flux density of 62.1 
mJy. The northeastern component has a peak flux density of 6.5 mJy beam~^ and an integrated flux density of 33.5 
mJy. Neither component has an optical counterpart in the deep optical data. 

J014439-002605 (CJ): A newly-identified variable source. The Stripe 82 data show two components which were 
previously unresolved in the FIRST survey. The central (variable) component has a peak flux density of 14.7 mJy 
beam~^ and an integrated flux density of 15.4 mJy. The component 4" to the east has a peak flux density of 0.7 mJy 
beam~^ and an integrated flux density of 2.0 mJy. The central source matches (within 0.03") to a stellar source of r 
— 21.1 in the deep optical data. The component to the east has no optical match. 

J020234+000301 (CL): Previously known to be variable, and also found to be variable in the Stripe 82 radio data. 
The three sources all have jet-like structures in the new Stripe 82 data. The central (variable) component has a peak 
flux density of 41.7 mJy beam~^ and an integrated flux density of 44.6 mJy. The lobe 12" to the north has a peak 
flux density of 4.6 mJy beam~^ and an integrated flux density of 6.6 mJy. The lobe 28" to the south has a peak 
flux density of 29.5 mJy beam~^ and an integrated flux density of 49.0 mJy. The central component coincides with a 
z = 0.366 quasar. The southern component is 2.95" from an r — 22.9 galaxy in the deep optical data, but since the 
morphology is very suggestive of a double-lobed radio galaxy, this is likely a coincidence. 

J021202-002749 (CL): Previously known to be variable. Three distinct components were already detected in the 
FIRST data, but the eastern lobe has now been resolved into two components of peak flux density 1.3/4.2 mJy beam~^ 
and integrated flux density 6.2/4.5 mJy closer/further from the center, respectively. The corresponding FIRST lobe 
had a peak flux density of 6.0 mJy beam~^ and an integrated flux density of 9.2 mJy. The western lobe has a Stripe 
82 peak flux density of 21.4 mJy beam~^ and an integrated flux density of 30.2 mJy, while the FIRST lobe had a peak 
flux density of 30.2 niJy beam^^ and an integrated flux density of 31.8 mJy. The deep optical imaging shows a source 
0.14" from the central (variable) component. It is classifled as a galaxy of magnitude r — 23.5. The other components 
have no optical matches. 

J021301-001814 (CJ): Previously known to be variable. Called a PS bv lde Vries et all ()2004[ ). the new data begin 
to resolve this source, showing evidence of CJ morphology. An optical source 0.13" from the variable radio source is 
classified as a galaxy of r = 22.2. 

J021553+001826 (CJ): Previously known to be variable. The source showed PS morphology in the FIRST survey, 
but the Stripe 82 data show evidence of a 2a extension to the south-west. The closest optical match (at 0.1") is a 
galaxy of magnitude r = 19.6 in the deep optical catalog. 

J021756-000936 (CJ): Previously known to be variable. This source was previously one component and classifled 
as CJ, but it has been resolved into two sources in the Stripe 82 data. The central (variable) component has a peak 
flux density of 2.4 mJy beam~^ and an integrated flux density of 6.7 mJy. The southeast component has a peak flux 
density of 0.5 mJy beam~^ and an integrated flux density of 1.4 mJy. A faint jet-like structure extends up to the 
northwest of the core. There appears to be an object approximately 1" from the core component in the deep optical 
catalog; however, an artifact in the image resulted in its not appearing in the catalog. 

J02 1840-001516 (PS): Previously known to be variable. This source is included in the list because Ide Vries et al.l 
([2004) classified it as a complex (CX) source, with three different components in an unusual configuration. They were 
unclear if the various components related at all. This multi-component source has been resolved into four separate 
components with the Stripe 82 data, and in the new data, the sources to the North look unrelated to the variable 
source. The component 14.8" to the northwest coincides (within 0.61") with a r = 23.7 galaxy in the deep optical 
data. There appears to be a double-lobed radio galaxy associated with this source. The variable radio source, which 
corresponds to a z = 1.171 quasar, is therefore a simple PS. 

J021945+001943 (CJ): A newly-identified variable source. The source consists of a central component with a 
jet-like feature extending to the south-west. The radio source matches to an optical quasar at z = 1.266. 

J220513-000428 (CL): A newly-identified variable source, and a clear double-lobed radio source. Three compo- 
nents (the core and two lobes) are already visible in the lower-resolution FIRST data. The Stripe 82 data resolves the 
source even further, with the southern lobe listed as two extended components with integrated flux densities of 118.2 
and 49.4 mJy. The northern lobe is extended as well, with an integrated flux density of 43.7 mJy, and the (variable) 
core is point-like with a peak flux density of 11.5 mJy beam~^. Within 0.1" of the core is a r = 21.8 object classified 
as a galaxy in the deep optical catalog. 

J220755-000215 (CL): Previously known to be variable. The FIRST catalog shows the central source and an 
extremely faint point source 21" to the southwest with a peak flux density of only 0.8 mJy beam"^. This source is 
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likely to be a side-lobe, and does not appear in the Stripe 82 image or catalog. In the Stripe 82 data, the central 
FIRST source has been resolved into two components (a central source and a lobe to the northeast with a peak flux 
density of 1.8 mJy beam~^ and an integrated flux density of 3.8 mJy), and the image implies the existence of a third 
component to the southwest that did not make it into the catalog. The source therefore appears to be a double-lobed 
radio source. The central (variable) component has an optical match 0.14" away in the deep optical catalog that is 
classified as stellar with a magnitude r = 21.1. 

J221909+003112 (CL): Previously known to be variable. Whereas FIRST saw only one component, we see that it 
has now been resolved into five different components in the Stripe 82 data: the component to the south of the central 
component has a peak flux density of 0.4 mJy beam^^ and an integrated flux density of 0.5 mJy; a component 2.4" to 
the north has a peak flux density of 0.7 mJy beam^^ and an integrated flux density of 0.8 mJy; a component 4.7" to 
the north has a peak flux density of 1.4 mJy beam^^ and an integrated flux density of 2.5 mJy; and the point source 
component to the northwest has a peak flux density of 0.4 mJy beam^^. It is unclear whether this last component is 
related. The central (variable) component matches (within 0.11") to a r = 22.4 object classified as stellar in the deep 
optical catalog. None of the other components have close matches in the deep optical catalog; 

J222235+001536 (CJ): Previously known to be variable. The source was classified bv ' de Vries et all (J20041 ) as 
having a CJ morphology, but we would have called this source a PS based only on the FIRST data. The Stripe 82 
data show more structure, giving more credence to a CJ classification. An optical match at 0.13" from the source is 
classified as a QSO at z = 1.362. 

J222729+000522 (CX): Previously known to be variable. The FIRST data show three components - a central 
source and two lobes - but the southwestern lobe has now been resolved into two separate components. In the FIRST 
data, this lobe had a peak flux density of 17.6 mJy beam~^ and an integrated flux density of 34.5 mJy. In the Stripe 
82 data, the southernmost component has a peak flux density of 6.1 mJy beam~^ and an integrated flux density of 
9.1 mJy. The northern component of that lobe has a peak flux density of 3.5 beam~^ niJy and an integrated flux 
density of 24.2 mJy. The north-eastern lobe has a Stripe 82 peak flux density of 6.1 mJy beam~^ and an integrated 
flux density of 46.3 mJy. The FIRST data gave a peak flux density of 28.8 mJy beam~^ and an integrated flux density 
of 45.5 mJy for that lobe. The deep optical data shows a source 0.22" from the central (variable) component. It is 
classified as a quasar with z — 1.518. None of the other components show close matches. A greyscale image of this 
source is shown in Figure [T3] to give another view of the complex morphology. 

J224224+005513 (CL): Previously known to be variable. The FIRST source has been resolved into two components 
in the high-resolution Stripe 82 data, with a third component (to the northwest) just below the detection threshold. 
The component to the south has a peak flux density of 1.2 mJy beam~^ and an integrated flux density of 2.3 mJy. 
The deep optical data show a source 0.17" from the central (variable) component that is classified as stellar with a 
magnitude r = 21.0. 

J230002+000046 (CL): A newly-identified variable source. The variable component corresponds to the core of a 
double-lobed radio galaxy. Neither of the lobes satisfy the variability criteria. The northwestern lobe had a FIRST 
peak flux density of 1.2 mJy beam"-'^ and an integrated flux density of 2.5 mJy, and in the Stripe 82 catalog it has a 
peak flux density of 0.4 mJy beam"-'^ and an integrated flux density of 3.0 mJy. The southeastern lobe had a FIRST 
peak flux density of 3.8 mJy beam^^ and an integrated flux density of 4.4 mJy, and in the Stripe 82 catalog it has a 
peak flux density of 2.2 mJy beam^^ and an integrated flux density of 5.1 mJy. None of the components have optical 
matches within 3". 

J231541+002936 (CJ): Previously known to be variable, and also found to be variable in the Stripe 82 data. 
This source had a point source morphology in the FIRST data, but has now been resolved into three components. 
The central component has a peak flux density of 17.0 mJy beam~^ and an integrated flux density of 23.8 mJy. An 
extended component 2.2" to the east has a peak flux density of 0.5 mJy beam~^ and an integrated flux density of 8.0 
mJy, and a component 7.0" to the southeast has a peak flux density of 0.4 mJy beam~^ and an integrated flux density 
of 1.7 mJy. The two point source features to the south and northeast are too faint to make it into the catalog, which 
has a flux threshold which depends on radial distance from bright sources. The central (variable) component matches 
to an optical quasar at z = 1.358. 
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TABLE 1 Epoch III Variable Sources 
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PS 


PS 
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8.31 


134 
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59 


0.35 


1.36 


PS 


PS 


6.60 


6.04 
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7.34 


7.13 


165 


6.74 


8.92 


126 


1.36 


1.35 


PS 


PS 


1.97 


1.24 


141 


1.95 


2.22 


148 


2.25 


2.88 


55 


1.02 


1.48 


PS 


PS 


1.41 


2.25 


134 


1.74 


1.61 


136 


1.90 


2.44 


51 


1.27 


1.73 


CL 


CL 


4.10 


3.95 


136 


3.68 


3.10 


136 


4.01 


5.24 


54 


1.23 
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PS 


PS 


1.73 


1.23 


139 


2.44 
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150 


2.23 


2.82 


54 


0.95 
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PS 


PS 



00 41 52.484 +00 03 02.16 
00 42 02.303 -00 16 06.35 
00 43 37.907 -00 12 44.20 
00 46 22.696 -00 10 22.33 
00 47 08.920 +00 06 03.02 
00 52 06.443 +00 14 10.16 

00 57 22.306 +00 16 44.67 

01 00 29.413 -00 24 24.05 
01 00 33.502 +00 22 00.17 
01 01 48.758 -00 06 46.92 
01 01 51.209 -00 02 32.85 
01 01 55.997 +00 41 29.20 
01 05 59.953 +00 19 04.29 
01 08 28.579 +00 03 12.57 

*01 11 06.789 +00 08 46.51 
01 12 48.603 -00 17 24.72 
01 14 34.309 -00 25 59.70 
01 21 57.479 -00 12 27.02 
01 22 47.743 -00 02 56.83 
01 29 56.720 +00 23 38.38 
01 31 00.177 +00 27 47.11 
01 31 51.644 +00 07 41.29 
01 32 20.294 -00 00 54.17 
01 33 16.865 +00 02 17.15 
01 34 00.142 +00 09 31.37 
01 44 39.861 -00 26 05.31 
01 52 14.399 +00 15 02.01 

01 59 40.965 +00 21 37.18 

02 00 06.382 -00 10 48.18 
02 01 26.777 +00 01 45.54 
02 01 55.957 +00 32 13.90 

*02 02 34.322 +00 03 01.83 
02 03 15.893 -00 14 32.54 
02 14 39.295 -00 24 05.36 
02 14 54.909 -00 21 00.89 
02 19 45.435 +00 19 43.31 
22 04 55.826 -00 01 47.25 
22 05 13.524 -00 04 27.77 
22 06 00.185 +00 22 15.44 
22 06 13.125 +00 42 24.45 
22 08 15.086 +00 42 53.60 

*22 08 22.892 +00 23 53.05 
22 09 22.864 +00 28 45.23 
22 16 24.950 -00 10 10.22 

*22 20 36.320 +00 33 34.17 
22 23 40.593 +00 01 37.16 
22 25 30.691 -00 10 31.01 
22 25 45.603 +00 40 37.78 

*22 27 04.246 +00 45 17.54 
22 27 23.275 -00 25 35.82 

*22 27 26.543 +00 10 59.25 
22 29 09.760 +00 04 44.03 
22 30 19.198 +00 50 46.84 

*22 30 47.467 +00 27 56.61 
22 33 17.849 +00 34 38.87 
22 33 24.886 +00 09 33.30 
22 37 24.664 -00 15 25.53 
22 37 30.754 +00 31 07.75 
22 40 39.443 +00 21 04.30 
22 44 32.256 -00 20 35.83 

*22 44 48.100 -00 06 19.65 
22 46 08.961 +00 32 14.43 
22 46 11.896 +00 32 32.82 
22 46 24.124 +00 46 24.49 

*22 46 27.685 -00 12 14.19 

*22 47 30.195 +00 00 06.42 
22 48 04.840 +00 32 52.42 

*22 49 22.295 +00 18 04.51 
22 55 37.623 -00 01 44.13 

22 59 33.637 -00 28 22.60 

23 00 02.369 +00 00 46.48 
23 00 08.229 -00 28 16.88 
23 01 34.862 -00 20 15.94 
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TABLE 1 - Continued 



Epoch I Epoch II Epoch III 
R,A Dec S^pk '5^int rms Spk S-mt rms -Spk -S^int rms Maj /var Morphology 
(mJy) (mJy) (^Jy) (mJy) (mJy) (^^Jy) (mJy) (mJy) (^^Jy) (") (FIRST) (S82) 

23 03 10.630 -00 01 34.97 2.81 2.61 137 2.64 2.44 136 3.04 3.73 50 0.98 1.41 PS PS 

*23 03 14.818 +00 00 52.34 1.72 1.29 136 3.85 3.36 135 3.12 3.69 51 0.78 2.24 PS PS 

23 04 23.074 -00 04 17.34 2.87 2.78 142 2.15 2.66 144 3.84 5.10 80 1.21 2.37 PS PS 

23 13 09.856 +00 08 05.58 2.71 2.26 137 3.46 3.28 140 2.88 3.90 50 1.08 1.44 PS PS 

23 13 17.125 -00 12 19.84 1.04 1.17 135 1.40 0.72 136 2.29 2.73 51 0.80 2.62 PS PS 

23 13 21.460 +00 38 12.34 2.87 2.82 135 3.73 3.78 133 3.51 3.84 57 0.58 1.34 PS PS 

23 13 22.737 -00 02 12.97 7.21 7.02 148 7.76 7.56 150 7.70 10.34 55 1.38 1.43 PS PS 

23 14 58.458 -00 22 31.80 1.19 1.15 136 1.92 1.65 134 1.91 2.52 56 1.07 2.12 PS PS 

23 15 39.664 -00 01 11.07 3.43 3.56 145 3.65 3.58 144 3.43 4.46 62 1.13 1.30 PS PS 

*23 15 41.657 +00 29 36.59 13.94 16.32 150 22.28 24.16 150 17.03 23.85 59 1.53 1.71 PS CI 

23 15 48.190 +00 07 21.33 6.05 6.32 148 5.99 6.00 147 5.54 7.80 64 1.52 1.30 PS PS 

*23 15 58.666 -00 12 05.44 6.27 6.09 141 3.65 4.08 141 4.21 4.98 54 0.99 1.72 PS PS 

23 16 07.791 +00 31 07.72 6.00 5.79 146 4.96 4.29 145 4.67 6.54 58 1.14 1.32 PS PS 

23 18 10.992+00 17 49.85 2.26 2.13 143 3.35 2.97 142 2.53 3.54 55 1.27 1.57 PS PS 

23 18 24.423 +00 30 51.88 3.82 3.76 142 2.81 2.95 141 3.26 3.84 54 0.83 1.37 PS PS 

23 20 11.604 +00 12 19.67 3.32 2.94 147 2.96 2.45 146 3.45 4.24 70 1.21 1.43 PS PS 
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TABLE 2 Epoch I/II Variable Sources 





Epoch I 




Epoch II 




Epoch III 












RA Dec 


•^pk 


■Spirit 


rms 


•^pk 


■Spirit 


rms 


•^pk 


^int 


rms 


Maj /„,, 


Morph 
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(mJy) 


(mJy) ( 


(MJy) 


(mJy) 


(mJy) 


(^Jy) 


(mJy) 


(mJy) 1 


(pJy) 


(") 


( 


;first 


) (S82) 


00 43 32.712 +00 24 59.84 


108.53 


111.57 


153 


122.43 


125.68 


151 


76.14 


104.13 


698 


1.44 


1.19 


PS 


PS 


00 48 19.124 +00 14 57.13 


89.40 


91.66 


150 


100.00 


101.17 


153 


74.97 


102.69 


549 


1.27 


1.12 


PS 


PS 


00 52 05.568 +00 35 38.11 


81.46 


87.95 


140 


35.07 


40.46 


142 


25.54 


32.51 


480 


1.15 


2.71 


CJ 


CL 


00 52 12.473 +00 09 45.22 


9.77 


10.45 


138 


12.15 
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TABLE 3 

Fractional Variability 



Fractional Variability This Work de Vries Galactic 
f N N N 

< 1.25 
1.25 - 1.50 
1.50 - 1.75 
1.75 - 2.00 
2.00 - 2.25 
2.25 - 2.50 
2.50 - 2.75 
2.75 - 3.00 
> 3.0 

Total 89 123 39 
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TABLE 4 
Optical Identifications 



SDSS Match Rate Galaxy Fraction Stellar Fraction 

Variable Sources 60% ± 5% 72% ± 6% 28% ± 6% 

Non-variable Control 44% ± 5% 92% ± 4% 8% ± 4% 



TABLE 5 
Quasar. Fraction 

Spectroscopic QSO Photometric QSO 



Number sources in area 2083 6986 

Variable Source QSO Fraction 10% ± 3% 13% ± 4% 

Control Sample QSO Fraction 3% ± 2% 2% ± 2% 



